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ABSTRACT 

Interacting binaries consist of a secondary star which fills or is very close to filling its Roche lobe, 
resulting in accretion onto the primary star, which is often, but not always, a compact object. In 
many cases, the primary star, secondary star, and the accretion disk can all be significant sources of 
luminosity. SIM Lite will only measure the photocenter of an astrometric target, and thus determining 
the true astrometric orbits of such systems will be difficult. We have modified the Eclipsing Light Curve 
code (jOrosz fc HauschUd t 2000) to allow us to model the flux-weighted reflex motions of interacting 
binaries, in a code we call reflux. This code gives us sufficient flexibility to investigate nearly 
every configuration of interacting binary. We find that SIM Lite will be able to determine astrometric 
orbits for all sufficiently bright interacting binaries where the primary or secondary star dominates the 
luminosity. For systems where there are multiple components that comprise the spectrum in the optical 
bandpass accessible to SIM Lite, we find it is possible to obtain absolute masses for both components, 
although multi-wavelength photometry will be required to disentangle the multiple components. In 
all cases, SIM Lite will at least yield accurate inclinations, and provide valuable information that will 
allow us to begin to understand the complex evolution of mass-transferring binaries. It is critical 
that SIM Lite maintains a multi-wavelength capability to allow for the proper deconvolution of the 
astrometric orbits in multi-component systems. 

Subject headings: astrometry — stars: binaries: close — stars: fundamental parameters — stars: 
individual (QZ Vul, Cyg X-l, SS Cyg, V592 Cas, Sco X-l, AR UMa) 



1. INTRODUCTION 

The field of stellar evolution has matured to the point 
that we generally understand the entire evolution of an 
isolated star if we know its initial mass. While many de- 
tails remain to be solved, it is unlikely that SIM Lite will 
provide a revolution in our understanding of the evolu- 
tionary processes for single stars of low and intermediate 
mass. The same cannot be said for binary star systems, 
where there are at present a number of outstanding is- 
sues that remain unsolved. Among the most difficult is 
how mass transfer proceeds during the common envelope 
stage that occurs in the post-main sequence lives of close 
binaries. The formation of nearly every mass transferring 
binary requires a common envelope phase where the sec- 
ondary star helps to strip the atmosphere of the primary 
while driving mass from the syste m, and shrinking the 
binar y's orbit. Current theory fc.f.. lPodsiadlowski et al.1 
2003) has a difficult time explaining the creation of sys- 
tems with massive black hole primaries (Mbh > 9 M Q ) 
and very low mass secondary stars (M2 < 0.7 M Q ). Sim- 
ilar problems exist across the mass spectrum of close and 
interacting binaries (IBs) , which for the purposes of the 
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current study include systems with a black hole, neutron 
star, or white dwarf (WD) primary and a non-d egenerate 
companion. For example, iHowell et al.l (|2001bT ) explores 
the well known 2-3 hr period gap observed in cataclysmic 
variable (CV) systems, which is postulated to be a re- 
sult of cessation of mass-transfer from the main-sequence, 
low-mass secondary star to t he white dwarf comp anion 
for periods between 2-3 hours. Howe ll et al.l (|2001bO runs 
binary population synthesis models that show this the- 
ory is only correct if the secondary stars in these systems 
are up to 50% oversized, and thus both components are 
really much less massive than has been typically assumed 
via standard main-sequence mass-radius-temperature re- 
lations. Thus, if accurate masses can be determined for 
these systems, their formation theories can be directly 
tested. 

Typically in binary star work, radial velocity curves 
are used to obtain masses for the components. While it 
can be quite simple to get a radial velocity curve for one 
component of an IB, (the secondary star if it dominates 
the luminosity, else typically the disk), it is difficult, or 
even impossible, to get such data for their black hole, 
neutron star, or WD primaries. Thus, along with the 
radial velocity amplitude of the secondary or disk, one 
must also assume a mass and know the binary orbital in- 
clination in order to solve for the component masses. It 
is especially difficult to determine the binary inclination. 
In a large number of IBs the secondary star can be quite 
prominent at infrared wavelengths, and since the object 
is distorted, it ex hibits ellipsoidal v ariations as it orbits 
the primary (c.f. iGelino et al.ll200l1 ). The amplitude of 
these variations are dependent on the orbital inclination, 



and to a lesser extent, the properties of the secondary 
star. However, even a small amount of contamination by 
the accretion disk introduces considerable uncertainty in 
this method, and nearly all interacting binary systems 
have some level of contamination from the accretion pro- 
cess. Thus, this technique of using infrared ellipsoidal 
variations to determine the inclination fails if one cannot 
ascertain the spectrum and level of the contaminating 
source. 

This paper will show how SIM Lite, with multi- 
wavelength, microarcsecond astrometry, will be the first 
mission capable of determining the multi- wavelength as- 
trometric orbits of IBs. The field of optical astrome- 
try has recently been making great progress from both 
ground-based observations, as well as is posed to make 
dramatic leaps with the upcoming space-based missions 
GAIA and SIM Lite. The limits of ground-based astrom- 
etry have been recently pushed by both CHARA and 
PRIMA/VLTI, but neither are capable of the kind of as- 
trometric measurements needed to probe IBs. CHARA 
has multi-wavelength capabilities , but can only provide 
angu lar resolution to ^200 /ias ([ten Brummelaar et al.l 
2005), which is far larger than the reflex motions to be 
discussed in this paper, and requires very bright tar- 
gets. PRIMA/VLTI will achieve ~30-40 /xas precision 
Ivan Belle et al.ll200"cl , which may be at the limit of use- 
fulness for these systems, but only in K-band, which does 
not help distinguish individual component orbits, as will 
be shown. The GAIA mission will provide astrometry 
for ^10 9 objects with 4 - 160 /ias accuracy, for stars with 
V=10-20 res pectively , and does posses multi- wavelength 
capabilities (jCaccia ri 2009). However, GAIA is a scan- 
ning satellite and cannot perform pointed or time-critical 
observations, and in fact only achieves this accuracy by 
averaging ~80 individual measurements, the individual 
errors of which range from 36 - L431 /ias, again for stars 
with V=10-20 respectively (|Mignarj2005l) . Thus, GAIA 
cannot provide the high-precision, time-critical pointed 
observations needed to study IBs. In contrast, SIM 
Lite will be able to point at any desired object for any 
length of ti me, providin g single-measurement accuracy 
of ~1 /ias (|NASAll200"l . As well, SIM Lite will have 
~80 spectral channels, spanning 450 to 900 nm, thus 
providing multi-wavelength, microarcsecond astrometry 
(NASA 2009). 

SIM Lite will offer the first chance at measuring as- 
trometric orbits for a large number of IBs, thus directly 
yielding inclinations and allowing for the precise mea- 
surement of the masses of both components in these sys- 
tems. Determining the inclination of a zero-eccentricity 
system from astrometry is as straightforward as deter- 
mining the ratio the semi-major and semi-minor axes of 
the projected ellipse on the sky. Ideally, one can derive 
accurate values for the masses of both c omponents from 
the a strometry of a single component (| Benedict et al. I 
12000ft . but to do so requires that one knows the orbital 
period, the semi-major axis of the apparent orbit, the 
parallax of the system, and the mass ratio of the two 
components. SIM Lite can provide the first three of 
these four parameters, but the mass ratios must be es- 
timated for many of the main systems of interest. How- 
ever, for binaries where there is a significant amount of 
light from more than one component, SIM Lite can de- 
termine absolute masses for both components. Since SIM 



Lite only measures the location of the photocenter of the 
system, the observed reflex motion of the system will be 
wavelength-dependent. For example, in a binary where 
one component is hotter than the other, the measured 
motion in the blue part of the spectrum will be differ- 
ent from that measured in the red, with shorter wave- 
lengths tracing the motion of the hotter component, and 
vice-versa. Thus, if one has at least two astrometric or- 
bits at different wavelengths, and one knows the ratio of 
luminosities in each bandpass, (i.e. the spectra of the 
components), then one could reconstruct the individual 
orbits of each component and obtain absolute masses for 
both. Since SIM Lite is currently designed to have ~80 
spectral channels, and since multi-color photometry and 
spectroscopy already exists for most IBs of interest, this 
technique should be quite feasible for most IBs. 

In the next section we describe the procedure and code 
used to model the reflex motions of IBs. In section 3 we 
examine the results for "proto-typical" systems, includ- 
ing both X-ray binaries and cataclysmic variables, and 
estimate the required observing time required by SIM 
Lite to obtain reasonably accurate parameters for each 
system. We also briefly investigate how the presence of 
disk temperature gradients and hotspots in the systems 
affect the derived orbits. In section 4 we present a full 
modeling of simulated SIM Lite data for one system, and 
show the precision and accuracy of recovered astrometric 
parameters and derived system masses. We summarize 
our results in section 5. 

2. THE MODELING PROCEDURE: THE REFLUX CODE 

refluxQ is a code that computes the flux- weighted as- 
trometric reflex motions of binary systems. At its core is 
the Eclipsing Light Curve (ELC) code, which is normally 
used to compute light curves of eclipsing binary systems 
(Orosz & Hauschildt 2000). Besides specifying the pri- 
mary and secondary stars, ELC allows for the inclusion 
of an accretion disk. As with other light curve model- 
ing programs, an array of physical effects are taken into 
account, such as non-spherical geometry, gravity bright- 
ening, limb darkening, mutual heating, and reflection ef- 
fects. The program can either use a blackbody formula 
for local intensities of the stellar components, or inter- 
polate from a large grid of NextGen model atmospheres 
(Hauschildt et al. 1999). ELC also allows for up to 
two hot or cool spots to be placed on each star, and on 
the accretion disk. Thus, ELC can reproduce nearly any 
binary system, including complicated systems such as 
cataclysmic variables, RS CVn systems, and Low Mass 
X-ray Binaries (LMXBs). Additionally, we have modi- 
fied the ELC code to allow for a mixture of blackbody 
and model atmosphere intensities for use among differ- 
ent system components, and to allow for the possibility 
of a free-free, or bremsstrahlung, accretion disk, which 

follows the form F\ oc T _ 5e _ iH'A" 2 , where T is the 
temperature in Kelvin, and F\ is the flux in power per 
unit area per unit wavelength, A. 

reflux takes input parameters for a specified system 
and feeds them to the ELC program, which generates 
an intensity map of the system at a specified phase and 

1 reflux can be run via a web in terface from 
http://astronomy.nmsu.edu/jlcough/reflux.html Additional 
details as to how to set-up a model arc presented there. 



wavelength, composed of N points evenly spaced on a 
grid around each star and the disk, if present, reflux 
then computes the system's center of light position using 
the formula 



A' 



{X,Y) = J2Fi'(?>v)i 



(1) 



i=0 



where (X, Y) is the system's center of light ( "photocen- 
ter"), with the center of mass located at (0,0), and Fj is 
the flux of a grid point i, located at (x,y). This is done 
for a complete orbit at 8 different wavelengths, currently 
chosen to be the standard UBVRIJHK bandpasses, 
with the astrometric reflex motion output in /zas using 
an estimate of the system's distance, reflux also cal- 
culates the observed spectral energy distribution (SED) 
of both the sum and the individual components (as seen 
at quadrature) for comparison to multi-wavelength pho- 
tometry to help constrain system parameters, reflux 
will also output a 3D animated gif of the system, the ap- 
parent multi- wavelength astrometric orbit, and the x and 
y components of this motion versus orbital phase. The 
actual light curves over an orbit are written to a text file 
for comparison with existing phased-resolved photomet- 
ric data. 

While SIM Lite only operates in the 4500 to 9000 
A bandpass, we include JHK photometry into the mod- 
eling process to help constrain the SEDs for systems of 
interest. We do this because it is often possible to detect 
the secondary star in the near-infrared, allowing one to 
better quantify its contribution at the wavelengths acces- 
sible to SIM Lite. (For consistency and possible compari- 
son to future work in NIR astrometry, we also include the 
J and K-band astrometric motions in the output plots.) 
Thus, before one can model the expected reflex motions 
for an IB, it is critical to have both reasonable param- 
eters for the system, as well as multi-color photometry. 
reflux requires as inputs estimates of the masses of the 
two stellar components, their temperatures, their radii, 
the orbital period, and orbital inclination. The eccentric- 
ity of the binary is assumed to be zero since in the ma- 
jority of these systems the companion star fills its Roche 
lobe, indicating a close orbit with strong tidal forces. In 
addition, the V magnitude and distance must be input 
for model normalization. If the system has an accretion 
disk, the temperature at its inner edge, the power-law 
index for the radial dependence of its temperature, and 
the radius of the inner and outer edges must be inputted. 
Alternatively, one can choose whether the disk follows a 
blackbody or free-free emission law. In the latter case, 
one has to specify and adjust a normalization constant 
for the luminosity of the disk emission to best match the 
observed SED. To demonstrate the use of reflux, we 
model several proto-typical IB systems below. 

3. MODELING THE REFLEX MOTIONS OF INTERACTING 
BINARIES 

It is critical that one uses reasonable system parame- 
ters to first match the observed SED before relying on the 
output reflex motions. There are three main scenarios 
for the visual SEDs of IBs: 1) systems where one stellar 
component dominates the optical SED, 2) where more 
than one component contributes to the optical SED, and 
3) disk-dominated systems. Besides the stellar compo- 



nents, and symmetric accretion disks, there are a num- 
ber of other features that are present in IB systems such 
as accretion disk hot spots, accretion streams, magnetic 
structures, jets, and other outflows that both affect the 
SED and which might have appreciable astrometric sig- 
natures. We investigate some of these features below, 
but to completely cover all of the behavior exhibited by 
IBs is beyond the scope of the current investigation. In 
the following we perform case studies for three different 
IB scenarios, using several well known objects. Addi- 
tionally, we examine how accretion disk or photospheric 
"hotspots" can distort the reflex motions of an IB. Fi- 
nally, for each system we also examine the feasibility of 
observing the system with SIM Lite, providing rough es- 
timates for the amount of time that will be required of 
SIM Lite to observe the system at a given precision, using 
the SIM Differential Astrometry Performance Estimator 
(DAPE) (|Plummerll2009]) . 

3.1. IBs with SEDs Dominated by the Primary or 
Secondary Star 

There are quite a number of IBs where the primary 
or secondary star completely dominates the SED of the 
system. For these systems, determining an astrometric 
orbit is straightforward. We examine two cases: 1) QZ 
Vul, a LMXB with a black hole primary and a cool, main 
sequence secondary star, and 2) Cyg X-l, a high mass X- 
ray binary (HMXB) with an O supergiant primary, and 
a black hole "secondary" . 

3.1.1. QZ Vul 

As shown in iGelino etHI (J2009I ). the SED of QZ Vul 
appears to be that of a reddened K2 dwarf from the 
optical through the mid-infrared. The system parame- 
ters are shown in Table Q] These parameters were in- 
put into reflux to model the astrometric motions, with 
a NEXTGEN model atmosphere used for the secondary 
star. The model SED, a 3D model of the system, and its 
wavelength-dependent reflex motions are shown in Fig.[TJ 
Since there is only a single visible component, QZ Vul 
does not exhibit any discernible wavelength dependency 
to its reflex motions. However, since the black hole is 
~24 times more massive than the secondary star, the 
K2 dwarf has a large apparent motion, and thus even 
at a distance of 2.29 kpc, the system's astrometric reflex 
motion (~8 /ias), in theory, would be detectable with 
microarcsecond astrometry. However, in practice, this 
particular system is so faint (V = 21.2), and has such 
a short period (0.3342 days), that, according to DAPE, 
SIM Lite can not reach the needed precision without in- 
tegrating for longer than the orbital period. Thus, one 
would need to find a much closer, and thus brighter, 
LMXB to observe with SIM Lite. The est imated orbital 
inclination for QZ Vul is 64° (|Gelinoll200il) , and the only 
effect one would see with a QZ Vul type system with a 
different inclination is a reduction or amplification of the 
y-component of the astrometric motion, corresponding 
to an increase or decrease of the inclination respectively. 

3.1.2. Cyg X-l 

Cyg X-l is a well-known HMXB dominated by an O- 
type supergiant orbited by a 10 M black hole with a pe- 
riod of 5.6 d. Cyg X-l exhibits a wide variety of behavior, 
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Figure 1. Top Left: SED plot for QZ Vul. Top Right: 3D model of QZ Vul at Phase=90°. Bottom Left: Reflex orbit for QZ Vul. Bottom 
Right: X and Y components of the reflex orbit versus phase for QZ Vul. 



Table 1 

Parameters for the QZ Vul 

System 



Parameter 




Value a 


Magnitude (V) 




21.2 


Distance (kpc) 




2.29 


Inclination (°) 




64.0 


Period (Days) 




0.3342 


Eccentricity 




0.0 


Mass of Star 1 


(Mq) 


0.32 


Mass of Star 2 


(M ) 


7.7 


Radius of Star 


1(R ) 


0.397 b 


Radius of Star 


2 (Rs) 


0.0 


T cff of Star 1 (K) 


4500 


T cff of Star 2 (K) 






a Values from Gelino (2001) 

b Secondary star fills Roche lobe 

such as highly variable X-ray and radio emission, some 
of which is presumably due to a relativistic jet. The op- 
tical light also varies, but much more weakly (Am « 0.05 
mag). Thus, in contrast to QZ Vul, the primary star will 
be the visible source in this system. We model this sys- 
tem using the parameters listed in Table[2j using a black- 
body for the O star. As can be seen in Fig. [21 there is no 
wavelength dependence to the astrometric reflex motion, 
and a single bandpass is sufficient to determine the or- 
bit. Even at a distance of 2.f kpc, the large separation of 
the components in Cyg X-l produces a significant astro- 
metric signature (~25 /zas), even though it is the more 



massive component that is responsible for the visible as- 
trometric motion in this system. Given the brightness of 
the system (V = 8.95), and the long orbital period, this 
system is easy to observe with SIM Lite. According to 
DAPE, 10 individual measurements, each with 2.5 //as 
precision, could be obtained in only a total of 1 hour and 
40 minutes of mission time, (given 5 minutes of target 
integration time per visit, broken into 5, 1- minute chops 
between the target and reference star.) This would pro- 
vide more than sufficient high-precision measurements to 
obtain a good solution for this system, directly yielding 
the inclination of the system and the semi-major axis of 
the primary star's orbit, and thus indirectly the masses 
of the components. For systems in which one component 
completely dominates the visual SED, having accurate 
photometry for the system is unnecessary in interpreting 
the astrometric data. 

3.2. IBs Where the SED is Not Fully Dominated by 
Either the Primary or Secondary Star 

3.2.1. SS Cygni 

SS Cyg is a bright, well-known, cataclysmic variable. 
In quiescence, SS Cyg has an apparent magnitude of V 
= 12.2. During outbursts, it brightens to V ~ 8.8. It 
has an orbital period of 6.603 h rs, and the secondar y 
star has a s p ectral type of K5 ()Harrison et al.l 120041) . 
iDubus et al.l (|2004f ) obtained near-simultaneous multi- 
wavelength photometry of SS Cyg. As shown in Fig- 
ure 5 of Harrison et al.l (|2007l) the white dwarf and ac- 
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Figure 2. Top Le/t: SED plot for Cyg X-l. Top Right: 3D model of Cyg X-l at Phase=90°. Bottom Left: Reflex orbit for Cyg X-l. 
Bottom Right: X and Y components of the reflex orbit versus phase for Cyg X-l. 

the luminosity of the system is completely dominated by 
the hot (10,000 K), optically thick accretion disk. 

We have generated a model system with the parame- 
ters listed in Table EH where the temperatures and radii 
of the stellar components have been derived from the 
literature. In the case of SS Cyg, we used a NEXTGEN 
model atmosphere for the secondary star, and assumed 
a blackbody for the white dwarf primary. We attempted 
both blackbody and free-free accretion disk models, and 
found that the free-free model provides the best match to 
the photometry. The results are shown in Fig. [31 where 
the model SED is compared to the photometry, with the 
contributions from each of the three components shown 
with separate symbols. As can be seen from these plots, 
the white dwarf dominates the flux contribution in the 
U-bemd, and fades thereafter, with the secondary dom- 
inating from B through K (though the accretion disk 
contribution becomes quite important in the infrared). 
A wire grid representation of the SS Cyg system at phase 
0.25 is also shown in Fig. [3j As shown in the bottom-left 
panel, the amplitudes of the reflex motions smoothly in- 
crease from B (semi-major axis a — 12 /ias) to R (a = 28 
/ias), and then decline at longer wavelengths due to the 
contribution of the accretion disk. The reflex motion in 
the t/-band is larger than in B, as it is almost completely 
dominated by the WD primary. In the -B-band, the two 
stars have similar luminosities, and thus the motions are 
a mixture of the individual astrometric orbits. If there 
was no accretion disk, and the white dwarf was invisible, 



Table 2 




Parameters for the Cyg X-l 


System 




Parameter 


Value a 


Magnitude (V) 


8.95 


Distance (kpc) 


2.10 


Inclination (°) 


35 


Period (Days) 


5.566 


Eccentricity 


0.0 


Mass of Star 1 (M Q ) 


17.8 


Mass of Star 2 (M ) 


10.1 


Radius of Star 1 (R ) 


16 


Radius of Star 2 (R ) 


0.0 


Tgg of Star 1 (K) 


32000 


T off of Star 2 (K) 





a Values from IHerrero et al.l 


( 1995) 





cretion disk dominate the blue end of the SED, while 
the secondary star becomes promi nent in the red and 
near-infrared. lHarrison et al.l (|2007|) modeled the SED of 
SS Cyg as a combination of the two stellar components 
plus a free-free accretion disk around the white dwarf. 
Due to the availability of simultaneous UBVRIJHK 
photometry, and well-constrained stellar parameters, SS 
Cyg is an ideal system to investigate what happens when 
more than one component in the binary system is visi- 
ble. As we demonstrate, SS Cyg has strongly wavelength- 
dependent astrometric reflex motions. We only consider 
the reflex motion in quiescence, since during outburst, 
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Top Left: SED plot for SS Cyg compared to photometry from IDubus et al.l (2001 ). Top Right: 3D model of SS Cyg at 
Bottom Left: Reflex orbit for SS Cyg. Bottom Right: X and Y components of the reflex orbit versus phase for SS Cyg. 



the total amplitude of the reflex motions for the sec- 
ondary star in this system would be a = 34.1 /xas. Thus, 
the true reflex motions are diluted, and one cannot deter- 
mine the actual astrometric orbit without modeling the 
system. However, with proper modeling, one can decon- 
volve the motion of the white dwarf + accretion disk and 
the secondary star, thus directly yielding the semi-major 
axis of each component's orbit, which when combined 
with the well-known period and distance to the system, 
directly yields absolute masses for each component, (see 
§4). 

Careful examination of the wavelength-dependent as- 
trometric orbits reveals that the actual center about 
which the motion occurs is wavelength dependent. Note 
that the C/-band reflex motion is centered on (0,0), but 
at the other wavelengths, the motions are offset from this 
position. This subtle effect has two causes: 1) the sec- 
ondary star in SS Cyg fills its Roche lobe, and therefore 
has a teardrop shape, and 2) gravity brightening for non- 
degenerate stars is significant. These two effects combine 
to offset the center of light of the secondary star from 
its center of mass. This creates an offset in the astro- 
metric motions in the bandpasses where the secondary 
star dominates the systemic luminosity. Since the WD 
is spherical, and does not have gravity brightening ef- 
fects, the C/-band motion remains centered at (0,0). For 
more detailed modeling on the wavelength dependence 
of the astrometric cen ter due to gravity brightening, see 
ICoughlin et all (l2010h . 



Table 3 

Parameters for the SS Cyg 



System 



Parameter 



Value a 



Magnitude (V) 

Distance (pc) 

Inclination (°) 

Period (Days) 

Eccentricity 

Mass of Star 1 (M Q ) 

Mass of Star 2 (M ) 

Radius of Star 1 (Rq) 

Radius of Star 2 (R ) 

T cff of Star 1 (K) 

T cff of Star 2 (K) 

Disk Inner Radius (Rq) 

Disk Outer Radius (Rq ) 

Disk Inner Temperature (K) 

Disk Temp Power-Law Exponent 



8.8 (max) 12.20 (min) 

159.5 

50.5 

0.275130 

0.0 

0.555 

0.812 

0.684 b 

0.015 

4400 

35000 

0.022 

0.407 

10000 c 

0.0 



a Values from lRobinson fc Bitiien | |2009D 

b Star fills its Roche Lobe 
c Disk is free-free 

One of the important consequences of these effects is 
that the actual orbital inclination is harder to derive. 
Naively, one would assume that the orbital inclination 
can be calculated from the ratio of the minor and major 
axes of the orbital ellipse. If we do this in the i?-band 
[i = cos _1 (17. 5/28.0)] we derive an orbital inclination 
of 51.4°, instead of the input value of 50.5°. Depend- 
ing on the bandpass, the orbital inclinations derived for 
our simulation of SS Cyg range from a minimum of 47° 



in the t/-band, to a maximum of 51.4° in the i?-band! 
While it is likely that these differences will be lost in the 
astrometric noise for most IB systems, there will be a 
small number of IBs (and many non-interacting binary 
systems) where this difference should be detectable, and 
proper modeling is required to accurately extract the sys- 
tem's inclination. 

With respect to observing this system with SIM Lite, 
according to DAPE, even in its low state with V = 12.1, 
a total mission time of 6 hours and 40 minutes would 
provide 10 measurements with respect to an astrometric 
standard star, each with 1.7 /ias precision and an indi- 
vidual integration time of 20 minutes per measurement, 
each composed of 20 chops between the target and ref- 
erence stars, with 1 minute exposures on the target and 
30 second exposures on the reference. As a minimum of 
7 points are needed to solve for a full astrometric orbit, 
(given that 7 parameters describe an orbit), and as the 
orbital period of SS Cyg is ~6.6 hours, this would suc- 
cessfully measure the full astrometric orbit with limited 
orbital smearing, (see Q. 

SS Cyg can be used as a template for cases of CVs 
and LMXBs where the secondary star is clearly visible 
in the optical. It simply requires changing the input val- 
ues of the stellar components, the relative prominence 
(and nature) of the accretion disk, and altering the or- 
bital inclination. Other changes, such as orbital period, 
mass ratio, and distance, act to scale the amplitude of 
the reflex motions. However, it remains critically impor- 
tant to have simultaneous multi-wavelength photometry, 
and a reasonable handle on the system properties to de- 
rive accurate astrometric orbits for these types of objects. 
As one would expect, the wavelength dependent behav- 
ior of the astrometric motions for these types of systems 
is more complex than for systems with only a single, 
visible component, but the potential scientific payoff is 
much higher. When varying the inclination angle for a 
SS Cyg type system, the y component of the astrometric 
motion decreases with increasing inclination angle for all 
wavelengths as expected, but at large inclination angles 
eclipse effects produce significant deviations from simple 
sinusoidal reflex motions. 

3.3. IBs with Accretion Disk Dominated SEDs 

The majority of LMXBs, and many CVs, have their 
optical light dominated by the accretion disk. The ac- 
cretion disk surrounds the compact object, and thus the 
observed astrometric motion is that of the compact ob- 
ject, although there are some complications arising from 
the extended structure of the disk. Here we examine two 
cases, the short period CV system V592 Cas, and the 
well-known LMXB Sco X-l. 

3.3.1. V592 Cas 

V592 Cas is a short-period (P or b = 2.76 hr) CV sys- 
tem consisting of a Roche lobe-filling M dwarf, a very 
hot white dwarf, and a disk that dominates the luminos- 
ity at optical and near-infrared wavelengths. The ma- 
jority of short period CVs closely resemble V592 Cas, 
and thus it can act as a prototype for disk-dominated 
CVs. The SEP of the system has been investigated by 
iHoard et al.l ([2009), who reproduced the observed pho- 
tometry using a model that included a red dwarf, white 



dwarf, and a two component blackbody disk around the 
white dwarf consisting of an inner flat component and 
an outer flared component. To model the mid- infrared 
excess, a circumbinary dust disk was also included. This 
dust disk is only important at mid-infrared wavelengths, 
and thus can be ignored in this present study. We model 
the reflex motions o f V592 Cas using the parameters from 
IHoard et al.l (|2009f ) shown in Table [4] while using a single 
component blackbody disk model with a flare inclination 
that is the average of the angles in their two component 
disk mod el. The SE P produced by reflux is compared 
to that of lHoard et al.l ((2009) in Fig. [H where the contri- 
butions from each component are also shown. We find a 
gr eater contribution o f flux from the M dwarf compared 
to IHoard et al.l (|2009f ). by a factor of ~l - 75. Th is is most 
likely due to the fact that IHoard et al.l ([2009) used the 
spectrum of a field M5.0V dwarf, while our model, which 
incorporates full Roche geometry, produces a star with 
the same effective temperature, but a larger surface area 
by the same factor of ~1.75, adopting the radius for an 
isolated M5.0V as 0.21 R su „, as shown by recent mea - 
surements of low-mass stars (c.f.. ILopez- Morales 2007). 
Either way, the M dwarf is a few percent of the total 
system flux, and thus does not greatly affect the derived 
astrometric measurements. 

As seen in FigureHJ due to the domination by the disk, 
the system shows less wavelength dependent motion than 
SS Cyg. However, there is still a discernible effect to- 
wards longer wavelengths as the secondary becomes a 
more significant source of flux. Thus, with enough mea- 
surements, one could in theory derive the individual com- 
ponent masses. The influence of the orbital inclination 
is readily apparent, and is uncontaminated by the sec- 
ondary star gravity effects that were visible in our models 
for SS Cyg. With respect to SIM Lite observing time, 
this system would require a moderate amount of SIM 
Lite mission time, given its small amplitude and short 
period. According to PAPE, if we constrain individual 
integration times to 15 minutes, (0.1 of the orbital pe- 
riod), to prevent orbital smearing, one could obtain 2.3 
fias precision per visit (composed of 20 target-reference 
chops of 45 second target integration each), which is on 
par with the amplitude of the observed reflex motions. 
Thus, in comparison to SS Cyg and Cyg X-l, discussed 
above, one might need more than 10 measurements to 
build up signal-to-noise and really nail down the system 
parameters. 30 measurements would take a total mission 
time of 17.5 hours, and thus we can say with reasonable 
confidence that this system would require less than a day 
of total mission time to observe so that the astrometric 
parameters could be obtained with reasonable precision. 

We decided to use this system to investigate the possi- 
bility of using the multi- wavelength astrometric capabil- 
ity of SIM Lite to discern two more subtle properties of 
disk-dominated systems: the disk temperature gradient 
and disk hotspots. The disk model employed in ELC al- 
lows for an inner disk temperature to be specified, along 
with a power law exponent, so that the disk temperature 
will vary radially according to the formula 



T(r) = T VI 



Ri 



(2) 



where T is the temperature at a given distance from 
the compact object, r, given a temperature Ti„„ er at 
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Figure 4. Top Left: SED plot for V592 Cas compared to models from lHoard et alj [f2009T ). Top Right: 3D model of V592 Cas at phase 
0.5. Bottom Left: Reflex orbit for V592 Cas. Bottom Right: X and Y components of the reflex orbit versus phase for V592 Cas. 



Table 4 

Parameters for the V592 Cas System 



Parameter 




Value a 


Magnitude (V) 




12.8 


Distance (pc) 




364.0 


Inclination (°) 




28.0 


Period (Days) 




0.115063 


Eccentricity 




0.0 


Mass of Star 1 (M ) 




0.210 


Mass of Star 2 (M ) 




0.751 


Radius of Star 1 (Rq) 




0.270 b 


Radius of Star 2 (R ) 




0.014 


T off of Star 1 (K) 




3030 


T off of Star 2 (K) 




45000 


Disk Inner Radius (Rq) 




0.0106 


Disk Outer Radius (R©) 




0.371 


Disk Inner Temperature 


(K) 


109,700 c 


Disk Temp Power-Law Exponent 


-0.75 


a Values from Hoard et al 


.1 (12009) 





b Star fills its Roche Lobe 
c Disk is blackbody 



the inn ermost radius, Rinner , an d a power law expo- 
nent x (|Orosz fc H auschildtJ 120001 ) . Usually a value of 
Y = — 0.75 is assumed for a steady-state disk (jPringld 
119811 ). but we wanted to test if changing \ to —0.5 or 
— 1.0, (which physically reflect a centrally irradia ted disk 
(|Friediung|[T98a IVrtilek et al.lll990t iBell fl999l) and its 
corresponding extrema), would have a discernible effect 
on the wavelength dependent astrometric reflex motion. 



Thus, we produced two models, one with x = —0.5 and 
the other with \ = —1-0, while also adjusting the inner 
disk temperature so that the resulting SED best matched 
that with x — —0.75. Neither the SEDs nor the reflex 
motions differed significantly from the models that used 
X = —0.75. Thus, it will not be possible to constrain 
the overall temperature gradient in accretion disks using 
astrometry. 

3.3.2. Sco X-l 

In many CV systems, the optical light is orbitally mod- 
ulated due to the presence of a hotspot on the outer edge 
of the accretion disk. This hotspot is where the accre- 
tion stream from the secondary impacts the disk. These 
spots typically have temperatures of 20,000 K, and can 
produce modulations on order of ±50% in visible band- 
passes (jMason et al.ll2002l ) . Because this feature is on the 
outer edge of the disk, it has the potential to distort the 
reflex motions in systems where it is prominent. In addi- 
tion, due to its higher temperature than the surrounding- 
disk, the hotspot could be the sou rce of strong line emis - 
sion, especially in HI or Hel lines (jSkidmore et al.i r2002) . 
Thus, its detection might be isolated, and the reflex mo- 
tions amplified, using very narrow bandpasses centered 
on the strongest emission lines. To investigate the ef- 
fects of a hotspot, we used the same V592 Cas model as 
above, but added a hotspot (20,000 K) on the outer edge 
of the accretion disk that leads the secondary by 30° in 
phase, and that is 30% of the U-band flux. The differ- 
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Figure 5. Top Left: SED plot for Sco X-l. Top Right: 3D model of Sco X-l at Phase=90°. Bottom Left: Reflex orbit for Sco X-l. 
Bottom Right: X and Y components of the reflex orbit versus phase for Sco X-l. 



Table 5 

Parameters for the Sco X-l System 



Parameter 



Value a 



Magnitude (V) 

Distance (kpc) 

Inclination (°) 

Period (Days) 

Eccentricity 

Mass of Star 1 (M ) 

Mass of Star 2 (M Q ) 

Radius of Star 1 (Rq) 

Radius of Star 2 (R ) 

T cff of Star 1 (K) 

T cff of Star 2 (K) 

Disk Inner Radius (Rq) 

Disk Outer Radius (Rq) 

Disk Inner Temperature (K) 

Disk Temp Power-Law Exponent 



11.1 (max) 14.1 (min) 

2.80 

50.0 

0.7875 

0.0 

0.42 

1.4 

1.473 b 

0.000 

4500 



0.026 

1.011 

100 000 c 

-0.75 



a Val ues f rom [S tccghs & Casarcs (2002), Gottlicb^TaD 
1T975T ), and Lasala & Thorstcnsen H 19851) 
b Star fills its Roche Lobe 
c Disk is blackbody 

ences between a model with a spot, and one without are 
extremely slight, with a maximum difference in models 
with and without a spot of 0.15 [i&s. Increasing the in- 
clination of the system does not have a significant effect. 
To be even remotely detectable, V592 Cas would have to 
be ~3 times closer, making the spot signature ~0.5 //as, 
although it would still only be ~5% of the reflex ampli- 
tude signature. Unless an IB has a hotspot that is much 



more dominant, such features will remain undetectable. 

Sco X-l is the proto-type LMXB consisting of a neu- 
tron star accreting from a low-mass companion. The 
X-ray luminosity of Sco X-l is very close to the Ed- 
dington limit for a 1.4 M neutron star. The sec- 
ondary star in Sco X-l has never been directly detected, 
though an orbital period of 0.787 days was detected 
through the an alysis of 85 years of visual photometry 
dGottlieb et al.lll975f). and from rad i al velocity measure- 
ments (jLasala fc Thorstensenlll985| ). Stccghs fc Casaresl 
(2002) detected narrow HI emission lines that they 
showed are from the irradiated secondary star. From 
these data they estimate masses for the components in 
this system of Mi = 1.4 M and M 2 = 0.42 M Q , thus 
making the secondary star a significantly evolved sub- 
giant. Since the accretion disk dominates the spectrum, 
we have simply assumed a Roche lobe-filling M type sec- 
ondary star corresponding to the observed mass, and an 
invisible neutron star. The stellar and accretion disk pa- 
rameters are listed in Table [5] and the accretion disk 
model we employ is similar to that for V592 Cas. 

As seen in Figure [5j the accretion disk dominates until 
the near-infrared when the secondary star begins to con- 
tribute. Thus, as also seen in Figure El the reflex motion 
has a weak wavelength dependence, but SIM Lite can- 
not observe in the near infrared where the reflex motion 
reverses, and thus it will be difficult to disentangle the 
components with such a weak wavelength dependence in 
the optical. However, it might be possible to recover the 
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Figure 6. Top Left: SED plot for AR UMa. Top Right: 3D model of AR UMa at Phase=90°. Bottom Left: Reflex orbit for AR UMa. 
Bottom Right: X and Y components of the reflex orbit versus phase for AR UMa. 



Table 6 

Parameters for the AR UMa System 



Parameter 




Value a 


Magnitude (V) 


14.50 (max) 18.00 (min) 


Distance (pc) 




85.0 


Inclination (°) 




65.0 


Period (Days) 




0.0805 


Eccentricity 




0.0 


Mass of Star 1 (Mq) 




0.70 


Mass of Star 2 (M Q ) 




1.40 


Radius of Star 1 (R Q ) b 




0.322 


Radius of Star 2 (R Q ) 




0.013 


T cff of Star 1 (K) 




3200 


T e fl of Star 2 (K) 




35000 


a Values from IHow 


ell et al.l 


(2001a) and 


IHarrison et al.l 112003) 





b Star fills its Roche Lobe 
c Disk is free-free 

contribution to the astrometric reflex motion for systems 
with slightly more prominent secondary stars (e.g., Cyg 
X-2). As seen in Figure [SJ Cyg X-l has a very small re- 
flex motion, only ~1.25 fias, and can get relatively faint, 
(V = 14.1 at its faintest), and thus this system is just at 
the limit of SIM Lite's capabilities. Adopting the faintest 
magnitude of 14.1 of the system to be conservative, and 
limiting ourselves to observations of 110 minutes, (0.1 of 
the orbital period), to prevent orbital smearing, we find 
via DAPE that we can achieve individual measurements 
with precisions of 2.2 /xas, almost twice the semi-major 



axis of the system. Thus, assuming we would need at 
least 50 measurements to begin to derive a reasonable 
orbit for this system, this system would require at least 
4 days of SIM Lite mission time, which may prove more 
expensive than the scientific payoff justifies. If we are 
able to observe in its high state when V = 11.1, we could 
achieve 30 measurements with 1.5 /xas precision in only 
20 hours of mission time, which is far more reasonable. 
When varying the inclination of a Sco X-l type system, 
the contribution from the secondary star becomes more 
significant as more of the disk becomes self-eclipsed and 
its total luminosity decreases, and thus alters the indi- 
vidual wavelength reflex motions accordingly. 

3.4. Polars and the Astrometric Signature of Magnetic 
Threading Regions 

3.4.1. AR UMa 

AR UMa is a "polar", a CV with a highly mag- 
netic WD primary an d a Roche lo be filling M5/ 6 dwar f 
([Harrison et al.ll2005l) . As shown in lSzkodv et at] (j 19991) . 
the U-band light curve is complex, being dominated 
by phase-dependent cyclotron emission. The near-IR 
light curve shows no rmal ellipsoidal variations, allowing 
IHowell e"t~ al. (2001a|) to estimate an orbital inclination 
of 70°. The magnetic field strength in AR UMa has 
been estimated to b e as high as B = 240 MG, though 
IHowell et all (|2001al) argue for a lower value of 190 MG. 
In these systems, at some point after material from the 
secondary star has passed through the inner Lagrangian 
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point, it is captured by the magnetic field of the white 
dwarf, producing a "magnetic threading region" that can 
be a significant source of luminosity in H Q . We investi- 
gate the possibility of astrometrically detecting this re- 
gion via the following modeling process. We model the 
system as a white dwarf + red dwarf system without an 
accretion disk, using the parameters shown in Table [51 
but for only the H Q bandpass add a narrow "accretion" 
disk halfway between the inner Lagrangian point and the 
WD. This disk then has a hot spot leading the inner La- 
grangian point by 45° in phase. The actual contribution 
from the disk has been made to be completely negligible, 
but the spot itself makes the system 25% brighter in H a . 
In effect, it is an isolated emission spot located between 
the WD and secondary star, mimicking a threading re- 
gion. 

As can be seen in Figure [6j there is a large amount of 
wavelength dependent astrometric reflex motion as ex- 
pected, ranging from ^10-20 /ias in the optical, with the 
WD dominating at short wavelengths and the secondary 
star dominating at long wavelengths. This should eas- 
ily allow the determination of both component masses 
with orbital coverage at a few passbands. According to 
DAPE, limiting individual observations to 12 minutes, 
(0.1 of the orbital period), to prevent orbital smearing, 
one could obtain individual measurements of ^6 /ias pre- 
cision at the system's brightest, (V = 14.5), or ^48 //as 
at its faintest, (V = 18.0). Thus, in its high state one 
could obtain 30, ^6 //as precision, measurements in as 
little as 9 hours of total mission time, but in its low 
state, assuming 100, ~48 //as precision, measurements 
would yield a meaningful solution, a total mission time 
of 30 hours would be required, which is still a reasonable 
amount of time for the scientific payoff in our opinion. 

Of particular interest is that in H a the threading region 
produces a detectable astrometric signature that has a 
phase offset of ^15° from the broadband wavelengths, 
due to the original 45° offset of the hot spot diluted by 
the light from the rest of the system. Thus, through 
use of narrow-band astrometry, one could recover the 
location of the threading accretion regions in these types 
of systems. When varying the inclination angle for a AR 
UMa type system, the y component of the astrometric 
motion decreases with increasing inclination angle for all 
wavelengths as expected, but at large inclination angles 
eclipse effects produce significant deviations from simple 
sinusoidal reflex motions. 

4. SIMULATED SIM OBSERVATIONS 

Although it is beyond the scope of this paper to 
perform a full simulation of exactly what parameters 
and to what precision one could extract with a given 
amount of SIM Lite time at specific wavelengths for 
each system, we can present a basic analysis for the 
SS Cyg system. To test the accuracy and precision to 
which we can recover system parameters, we generated 
simulated astrometric data for SS Cyg in the U and V 
bands, using the parameters listed in Table [3l but with 
an inclination of 54.0°, comprised of 10 measurements 
at random phases with 1.7 //as Gaussian noise added 
to each measurement. We then simultaneously solved 
for the astrometric parameters of period, semi-major 
axis, inclination, the position angle of the line of nodes, 
and the angle in the plane of the true orbit between 



the line of nodes and the major axis, while fixing the 
the eccentricity and the epoch of passage through 
periastron to be zer o, using the GAUSSFIT program 
(Peffervs et al.1 I1988D . with a procedure described in 
iBenedict et al.l (|2001| ). We recover inclinations of 55.6 
±5.1 and 54.0 ± 3.8 degrees, and semi- major axes of 
13.58 ± 0.94 and 23.58 ± 1.28 /xas, each for the U and 
V-bands respecti vely. Given the di stance to the system 
of 159.5 parsecs ()Dubus et al.H2004T ). and assuming that 
58% of the U-band astrometric motion results from the 
white dwarf, and 69% of the V-band flux results from 
the secondary star, as per the modeling, (see §3.2.1), 
this yields masses of 0.559 ± 0.070 and 0.815 ± 0.096 
M for the secondary star and white dwarf respectively, 
compared to the input masses of 0.555 and 0.812 
Mq. Thus, one will be able to disentangle individual 
component masses with reasonable uncertainties for this 
system with only ^8 hours of SIM mission time. 



5. CONCLUSION 

We have presented a modeling code, reflux, that 
is capable of modeling the multi-wavelength astromet- 
ric signature of interacting binaries. We have presented 
models for multiple IB systems, and find that for primary 
or secondary star dominated systems, the contamination 
of the photocenter is minimal and SIM Lite will provide 
good astrometric orbits if the system is bright enough, 
yielding absolute inclinations and indirect masses. For 
mixed component systems, SIM Lite will be able to di- 
rectly determine absolute masses for both components 
when multi-wavelength astrometric curves are obtained 
and combined with existing spectroscopy and multi-color 
photometry. For disk-dominated systems, SIM Lite will 
only be able to obtain astrometric orbits for a few sys- 
tems, due to the small relative motion of the compact 
object that the disk surrounds, but these should yield ac- 
curate inclinations. We find that while multi-wavelength 
SIM Lite data will be unable to distinguish between var- 
ious disk temperature gradients and disk hotspots, it 
should be able to determine the location of magnetic 
threading regions in polars, and thus similar effects in 
other systems, via narrow wavelength astrometry. In to- 
tal, SIM Lite should contribute greatly to our under- 
standing of interacting binary systems, especially if its 
multi- wavelength capability is maintained in the final de- 
sign. 
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